Multiple charging of InAs/GaAs quantum dots by electrons or holes: addition 

energies and ground-state configurations 
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Atomistic pseudopotential plus configuration interaction calculations of the energy needed to 
charge dots by either electrons or holes are described, and contrasted with the widely used, but 
highly simplified two-dimensional parabolic effective mass approximation (2D-EMA). Substantial 
discrepancies are found, especially for holes, regarding the stable electronic configuration and filling 
sequence which defies both Hund's rule and the Aufbau principle. 
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I. INTRODUCTION 

One of the most spectacular aspects of quantum-dot 
physics is that dots can be controUably charged by either 
electrons or holes and that one can measure, for each of 
the many charged-states, both the electronic spectrum 
and the charging energies. This is afforded either by in- 
jecting from a tip of a scanning tunneling microscope^ 
or by various gate structures i^*^*^i^iSiL2i2iiS4ii Since the 
energy scale of both single-particle levels and Coulomb 
interactions in quantum dots (QDs) (10^^ - 10^^ Ry) are 
a few order of magnitudes smaller than those of the real 
atoms ( ~ 1 Ry), dots can be loaded by as many as six* to 
teni^ electrons in colloidali^ and self-assembled^iSiS dots 
having confining dimension of < 50 A, and up to hun- 
dreds of electrons in larger 500 A electrostatically con- 
fined dotsjSj^ The "charging energy" /i(iV) is the energy 
needed to add a carrier to the dot that is already loaded 
by iV — 1 carriers, 



H{N) = E{N) - E{N - 1) 



(1) 



where E{N) is the correlated, many-body total energy of 
the ground state of the iV-particle dot. The "addition en- 
ergy" A (TV — 1, A*") (analogous to the difference between 
ionization potential and electron affinity) indicates how 
much more energy is needed to add the iVth carrier com- 
pared to the energy needed to add the {N — l)th carrier: 

A(iV-l,iV) = ijl{N) - ^i{N - I) 

= E{N) - 2E{N - 1) + E{N - 2) . (2) 

The typical electron addition energies for electrostatic 
dots^i^ are about 1-8 meV, and the stable spin- 
configuration follows the rules of atomic physics; that 
is, the s, p, d, ... shells are occupied in successive order 
with no holes left behind (Aufbau principle) and with 
maximum spin^ (Hund's rule). Recently, it became pos- 
sible to load and measure electrons"'-, and holes^tii^i^ 
into much smaller, epitaxially grown self-assembled dots 
of InGaAs/GaAs, where electron addition energies are 
about 10 - 60 meVr^ and hole addition energies are be- 
tween 10 - 30 meViii*i^ More interestingly, while elec- 
trons still follow the Aufbau principle, recent hole charg- 



ing experimentii show that holes have unusual charging 
patterns that defy the Aufbau principle and Hund's rule. 

Despite the importance of the dot charging problem 
and the great success achieved in experimentally record- 
ing the charging spectra, the theoretical understanding of 
charging and addition energies is still preliminary. Most 
theoretical works in the area were based on particle-in-a- 
potential modelji^ii^iiSiii neglecting interband (e.g., F- 
r) coupling, intervally (T-X-L) coupling effect, and the 
true atomistic symmetry (e.g., C2v for lens of zincblende 
material) which is lower than the shape symmetry. The 
most often used potential in such approaches is the 2- 
dimensional (2D) parabolic form, in which all of the 
above noted electronic structure effects are replaced by 
an effective mass approximation (EM A). In this 2D-EMA 
modelji^ the single particle levels have equal spacing, 
which equals a harmonic oscillator frequence u. Because 
of the simplicity of the model, all the Coulomb inte- 
grals can be related analyticalljii^ to a single s orbital 
Coulomb energy Jss, and therefore the addition energies 
are determined entirely by w and Jss ■ Although the 2D- 
EMA model can be attractive because of its algebraic 
simplicity, and availability of fitting parameters, actual 
self- assembled dots are significantly different from the 
description of EMA model, manifesting inter-band cou- 
pling and inter-valley coupling, strain effects, low atom- 
istic symmetry, as well as specific band offset profiles, all 
neglected by the 2D-EMA. It indeed has been recently 
shown^^ that hole ground state configurations predicted 
by the 2D-EMA model are qualitatively different from 
those measured by hole charging experimentsiii*i^ 

An atomistic pseudopotential description of electronic 
structure effects can be used instead of 2D-EMA to 
calculate charging energiesiiSiiSiSSiSiiS^ Here, we show 
that such an atomistic theory correctly reproduces the 
many-particle configurations as well as addition spectra 
for carriers in self-assembled quantum dots. We study 
systematically the electronic structure of self-assembled 
InGaAs/GaAs quantum dots, provide detailed infor- 
mation on the electron and hole single-particle spec- 
trum, many-particle charging and addition spectrum, as 
well as ground state configurations using single-particle 
pseudo-potential and many-particle configuration inter- 
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action (CI) methods. 

The rest of the paper is arranged as follows. In SeclTTl 
we introduce the basic concepts of charging and addi- 
tion energies, and show how to calculate these quantities 
in the single-particle pseudopotential plus many-particle 
CI scheme. In Sec. IIIII we give detailed results calcu- 
lated from pseudopotential-CI scheme the single-particle 
levels, Coulomb integrals and the ground state configura- 
tions as well as the addition energies. We contrast these 
results with 2D-EMA model. We summarize in Sec. IIVI 



step is formulated as, 
1. 



V'i(r) = ei , (3) 



where T4xt(r) is the external ("bare") potential experi- 
enced by the electrons or holes, and t4cr(r) is the screen- 
ing response. The single-particle orbital {ipi} and en- 
ergies {ei} are used in the second step to construct the 
many-particle wave functions {5'} and energies {E} from, 



II. THEORY OF DOT CHARGING AND 
ADDITION ENERGIES 

A. General equation for dot charging in the 
configuration-interaction approach 

The calculation of the total-energy E{N) of A^-particle 
dot requires obtaining first the single-particle states from 
an effective Schrodinger equation, and then the many- 
particle state from a many-particle treatment. The first 



E{N) = {^n\H\-^n) (4) 
where, the many-body Hamiltonian is, 



ijkl CTl ,fT2 (T3 ,(T4 



and, 



p?Cri,JO-2 



E 



drdr- 



, IPl a^ (l": S2)V'fc^3 . 82)^';^, (r, Sj) 

e(r — r')|r — r'| 



(6) 



are the screened Coulomb and exchange integrals. In the 
above Eqs. (|SJ and (|HJ, "cr" is a pseudospin index, i.e., an 
index of Kramers degenerate states, while "s" is the in- 
trinsic electronic spin. For electrons in InAs/GaAs QDs, 
the spin-orbit interactions is extremely small and can be 
neglected. In this case, the pseudospin a and intrinsic 
electronic spin s are equivalent. However, for holes, which 
have a mixture of heavy-, (H) light-hole (LH) and split-off 
character, an eigenstate of a has both s =] and s =[ com- 
ponents. The A^-particle wave functions can be solved 
using e.g. configuration interaction (CI) method^S^ by 
expanding the A^-electron wave function in a set of Slater 
determinants, |'I>ei,e2,- - ) = V'liV'la •••V'Ll^o), where 
creates an electron in the state . The v-ih many- 
particle wave function is then the linear combination of 
the determinants, 

l*7V^) = E ^^(ei, 62, • • • , bat) |$ei,e2,- ,e„) ■ 

(7) 

Once 1 4* at) is known, we can then calculate the corre- 
sponding total energies for the ground states as well as 
excited states using Eq. @. Once we solve the CI prob- 
lem, we get the order of total CI energy for various holes 
or electron configurations, so we can see if Hund's rule or 
the Aufbau principle or spin-blockade occurs. For exam- 
ple, Hund 's rule states that degenerate single-particle lev- 



els are occupied with maximum number of unpaired elec- 
trons, while the Aufbau principle states, non-degenerate 
single-particle levels are occupied in order of increasing 
single-particle energy. 

We construct all possible Slater determinants corre- 
sponding to N electrons or TV holes [i.e., we ignore the ex- 
citonic (electron-|-hole) excitations], using only the bound 
states of the dots, (i.e., we neglect all continuum states). 
The underlying electrons that are not considered explic- 
itly by this approach are represented by the dielectric 
screening function e(r — r') in Eq. ©. 



B. The Hartree-Fock equations for charging and 
addition energies 

The addition energies at CI level can be written as 
the Hartree-Fock (HF) addition energies plus the corre- 
lations, i.e., 

Aci(iV-l,iV) - AHF(A^-l,iV) + Aeorr.(^-l,^) , (») 

where Acorr. is the correlation energy correction to the 
addition energy calculated in HF. Since the HF equations 
are used by many experimentalists to deduce Coulomb 
energies f^iSiii we review it here. In the Hartree-Fock ap- 
proximation, where the effect of correlations is neglected 
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but the direct Coulomb and exchange interactions are re- 
tained, simple expressions can be derived for the addition 
energies. The total energy of N electrons is simply, 

occ. occ. 
EkF = ^^i(T + ^ {Jiaja' — Ki„j„i) , (9) 

where, i is the single-particle level index of all occupied 
states, and e.^ is the corresponding single-particle energy. 

The J,a,3a' = rp;;^ and K^^,^^' = in Eq. © are 

Coulomb and exchange integrals, respectively. Since the 

I 



From above equations, we see that to predict A(7V— 1, TV) 
for the first two cases (A^=3), we need to know 4 param- 
eters (Jss, Jspi, Ksp^, - e^), but for A^=6, we need 
to know 9 parameters. Because of the large number of 
parameters needed, the analysis of charging effects in the 
literature resort to additional approximations aimed at 
reducing the number of parameters using simplified ef- 
fective mass models. Equations (|12|) are going to be used 
below to contrast the charging spectra deduced from sim- 
plified literature models vs. our more complete treat- 
ment. 



spin index "cr" is not an actual electronic spin, but rather 
an index for two Kramers degenerate states, in principle 
the exchange integrals Kia-ju' are not simply diagonal 
in tJ, a', as has been widely used in the dot charging 
literature i^iiitiSiiSi However, adopting the literature ap- 
proximation, Kia-ja'=Ki,j baa' and considering the s and 
two p orbitals, ep^ < t-p^ and assuming the particle filling 
order follows Hund's rule, as shown in FigQ the total 
energies for iV=l, 2, 3, 4 electrons in the Hartree-Fock 
approximation are 



(10) 



(11) 



(12) 



C. Atomistic treatment of the single- particle 
problem 

The most general treatment of the single-particle prob- 
lem of Eq.© describes both 14xt and V^cr atomistically, 
much in the same way as molecules are treated quantum- 
mechanically. In this description Vc^t{v) is a superposi- 
tion of the ionic potential of individual atoms of type a 
at lattice site n, 

^ext(r) = E E ^io"n (r - - Rn) , (13) 



£;hf(i) = cs , 

£;hf(2) - 2e, + J,,, , 

-Ehf(3) = leg + Cpi + Js,s + 2Js^pi — '^s,pi , 

-^Hf(4) — 265 -\- -t- ~1~ J S.S ~t- s.p\ ~t~ 2(/s^p2 ~t~ '^Pl,P2 ^S,pi ^S,P2 ,P2 

Ehf{5) = 2es -I- 2ep^ -1- Cp^ + Jg^s + 4Js,pi -I- ^Jg^p^ + Jpi,pi + 2Jp^^p^ — IKg^p^ — K, 

-^Hf(6) — 265 -\- 26p^ -\- 26p2 ^ Js,S ^ ^Js,pi ^'^S.p2 '-^Pl.pi ^ '-^P2,P2 ~^ ^*^Pl,P2 

—2K —IK —IK 

'^J-^S.pi '^-^^S,P2 ^-^^Pl.P2 ' 

We can then readily calculate the charging energies Eq. Q in this approximation. 
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Similarly, we can calculate the addition energies Eq. Q as follows, 
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FIG. 1: The sequence of electron configurations for filling 
A'^ electrons to the quantum dots according to both 2D-EMA 
model and pseudopotential calculations. Note that the p lev- 
els are split (into pi and P2), so are the d levels (di, d2, ds). 

and v[^^^ is —Z^jr in an all-electron (core-t- valence) treat- 
ment (where Za is the atomic number) or v'^ iv) in the 
pseudopotential (valence-only) scheme (where (r) is 
the ionic pseudopotential. The potential 14xt naturally 
contains the correct point-group symmetry of the object, 
through the atomic position vectors {tc(,R„}, and in- 
cludes atomic relaxation if appropriate (again, through 
the atomic positions), as well as chemical inhomogene- 
ity (alloying) or surface-passivation effect. The screening 
response 14cr(r) of Eq. Q is in general a functional of 
the density-matrix p{y, r') and can be described e.g. via 
Hartree-Fock or the density functional theory, both re- 
quiring a self-consistent (iterative) solution to Eq. Q. 
These approaches are currently limited to small dots, 
relative to the 10^ - 10^ atom dots which charging ex- 
periments exist. Furthermore, LDA suffer from the fa- 
mous "LDA error" , whereby the band gap and effective 
masses are badly underestimated. Higher-order method 
such as GW approximation or time-depend DFT have 
yet to demonstrate applicability to large dots for which 
high-quality experiments exit. 

An approximation to the screening V^cr[p(r, r')], which 
allows calculation on large dots, and fixes the LDA error, 
is provided by the "screened pseudopotential approach" , 
where it is assumed that 14cr(r) can be described as a 
superposition of screening potentials wicr (r) of the in- 
dividual atoms, and lumping together v\^^(v) -I- wi"r (r) 

to yield a screened atomic pseudopotential wlpmlr), such 
that, 

K=xt(r) + T/scr(r) = F,o(r) - - Rn) . 

Q n 

(14) 

Here, i'epm(i") is determined semi-empirically. Unlike the 
classic empirical pseudopotential methodf24 which fitted 
only to eigenvalues, here we require that when Eq. (PJ is 
applied to the underlying hulk periodic solids containing 
atom {a}, {wcpm} reproduce the measured band ener- 
gies, effective-mass tensors, deformation potentials, and 
the single-particle wave functions have a large overlap 
with the corresponding LDA wave functions^--'— In Eq. 



(|14|l . a nonlocal potential V^o(r) is also added to the total 
potential to represent the spin-orbit interaction. In our 
approach, the potential of an As atom depends on the 
number of Ga and In atoms around it as 

71 4 Ti 

?;As(Ga„In4_„) = -WAs(Ga) -\ —vas{^'o) , (15) 

where n is the number of Ga atoms around the As atom. 
In this atomistic approach, one assume that Wcpm(r) is 
transferable to different environments. Note that a fixed 
Va is a good approximation if the dot has no free sur- 
faces (as is the case in self-assembled dots, where only a 
strained interface between chemically-similar materials is 
present). For surface atoms in free-standing dots, a sep- 
arate t;opm(r) is fittedSi to LDA surface calculations. For 
InAs/GaAs dots, we use the pseudopotentials of Refill 
These pseudopotentials have been tested not only for the 
InAs and GaAs binaries, but also for alloys and super- 
lattices of the corresponding ternaries.-^ 

Once {wopm(i")} is known, one can solve Eq. ||2Jl for 
the bulk solid, quantum wells superlattices, quantum- 
wires or quantum dots by adopting a supercell approach 
where the respective objects are placed. In our case, 
Eq. ||3J) is solved using the "linear combination of Bloch 
bands" (LCBB) method^^ where the wave functions i>i 
are expanded as, 

^^w^EECCeW. (16) 

n,k A 

In the above equation, {0|^^],„(r)} are the bulk Bloch 
orbitals of band index n and wave vector k of ma- 
terial A (= InAs, GaAs), strained uniformly to strain 
e*. The inclusion of stain-dependent basis functions im- 
proves their variational flexibility. We use e = for the 
(unstrained) GaAs matrix material, and an average e 
value from VFF for the strained dot material (InAs). For 
the InAs/GaAs system, we use n = 8 (including spin) for 
electron states on a 6x6x16 k-mesh. Note that the po- 
tential V^xt (r) + Ksrc (r) contain full strain effects through 
the use of relaxed atomic positions, in addition to the 
explicit strainiS and alloy composition'^^ dependence. 

In the atomistic approach to the single-particle prob- 
lem, one includes (i) multi-band coupling [different n in 
Eq. (4)]; (ii) inter- valley coupling [different bulk fc-points 
in Eq. (4)]; (iii) spin-orbit coupling \Vso in Eq. ((1^ ]: 
(iv) the proper strain profile [by relaxing {tq,, R„} in Eq. 
(|14|l to minimize strain] ; (v) realistic chemical profile [dis- 
tributing the species {a} as in a random alloy, or inter- 
diffused interfaces'^] . The ensuing single-particle orbitals 
{ipi} transform like the representations of the point-group 
created by the ionic positions. These underlying atom- 
istic structures could break the symmetry represented 
by the macroscopic shape of the quantum dots. For ex- 
ample, a lens-shaped dot has a macroscopic cylindrical 
symmetry with [IIO] and [110] being equivalent, but if 
the dot is made of a zincblende material, the real sym- 
metry is C2v^ where, [110] and [110] are not equivalent. 
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Yet the continuum models do not "see" the atomistic 
symmetry (the "farsightedness effect"—). Therefore, as 
discussed in Ref. |^ in reahty the atomistic wave func- 
tion need not to be simple "pure" s-like or p-like. As 
a result, the Coulomb energy J^j and exchange energy 
Kij obtained with atomistic single-particle orbital {"tjji} 
do not have simple relationships^^ as predicted by con- 
tinuum 2D-EMA model. The deviations of the atomistic 
calculated Js from the simplified 2D-EMA ones are going 
to lead to new physical behavior (e.g., new ground state 
symmetry of the many-particle state), as illustrated be- 
low. 



D. Continuum treatment of the single-particle 
problem: EMA 

It is sometimes customarji^iiiiiSiS^ to avoid an atom- 
istic description of the single-particle problem in favor 
of a single-band particle-in-a-box model. In this ap- 
proach, one sets 14cr=0 and replaces 14xt of Eq. (3) 
by a pure external potential, describing the macroscopic 
shape of the object, e.g., a box represents a quantum 
well, or a sphere with finite or infinite barriers repre- 
sents a quantum dotii^iiSiS In some cases, this V^xt is 
calculated "realistically" from a combination of band off- 
set, the gate potential and the ionized impuritiesr^i^^ 
but it is treated nevertheless as a macroscopic external 
field. Under this approximation, simple results can be 
obtained for cylindrical QDs, where the angular momen- 
tum L and are good quantum numbers. The elec- 
tron (hole) single-particle levels have well defined shell 
structures, with non-degenerate s shell, two-fold degen- 
erate p shell, and three-fold degenerate d shell, etc. Fur- 
thermore, for parabolic confining potentials, the single- 
particle levels have equal spacing between two adjacent 
shells, e.g. ep — es=ed — ^p- Another reason for the attrac- 
tiveness of the single-band particle-in-a-box approach to 
the single-particle problem [Eq.©], is that the many- 
particle problem [Eq.Q] becomes simple. For example, 
the quantum Monte Carlo (QMC) approach is currently 
applied only within the single-band particle-in-a-box ap- 
proach for such large objects as self-assembled quantum 
dots.^^ In a 2D-EMA model all the Coulomb integrals 
needed for charging and addition energy can be all re- 
lated tc^^ Jss, 

J - ^■''^ J 

K^p = c%\jss- (17) 

For example, Jgp = 3/4Jss, Ksp = l/4Jss. Therefore the 
charging/addition energies and ground state configura- 
tions are totally determined by Jss and the single-particle 
energy spacing oj. 

However, real self-assembled quantum dots grown 
via the Stranski-Krastanov techniques, are not well- 
described by the single-band particle-in-a-box ap- 
proaches, despite the great popularity of such approaches 



in the experimental literaturesiiiii^iiSiS The model con- 
tains significant quantitative errors^S and also qualitative 
errors, whereby cylindrically symmetric dots are deemed 
to have, by symmetry, no fine-structure splitting, no po- 
larization anisotropy, and no splitting of p levels and d 
levels, all being a manifestation of the "farsightedness 
effect"^ 



III. RESULTS 

Using the pseudopotential approach for single-particle 
and configuration interaction approach for the many- 
particle step, we studied the electron or hole addition 
energy spectrum up to 6 carriers in lens-shaped InAs 
dots embedded in a GaAs matrix. We study dots of 
three different base size, h— 20 nm, 25 nm and 27.5 nm, 
and for each base size, two heights, h=2.b nm and 3.5 
nm. To study the alloy effects, we also calculated the 
addition spectrum for alloy dots Ini_2:Ga2:As/GaAs of 
h/b— 3.5/25 nm dots, with Ga composition x=0, 0.15, 
0.3 and 0.5. In this section, we give detailed results of 
the single particle energy levels and Coulomb integrals, 
and the addition energy as well as ground state configu- 
rations. We also compare the results with what can be 
expected from 2D-EMA model. 

A. Single-particle level spacing: Atomistic vs. 
2D-EMA description 

1. Electron levels 

We depict in Fig. [5] the calculated energy-level dia- 
gram of a pure lens-shaped InAs/GaAs quantum dot, 
with height h=2.5 nm and base b=20 nm. Figure[21shows 
that the electron confinement energy is 230 mev, some- 
what larger than the hole confinement energy (190 meV). 
The p levels are split as are the d levels, even though the 
dot has macroscopic cylindrical symmetry (see below). 

The pseudopotential calculated electron single-particle 
energy spacings are summarized in Table ^ for QDs of 
different heights, bases, and alloy compositions. Tabled 
gives the fundamental exciton energy Ex calculated from 
CI approach for each dot. These exciton energies are 
between 980 meV - 1080 meV for pure InAs/GaAs dots, 
and can be as large as 1297 meV for Ini_a;Gaa;As/GaAs 
alloy dots. This range agrees very well with experimental 
results for these classes of dots, ranging from 990 meV to 
1300 meV,«iiii38,39 

(a) s-p and p-d energy .spacing: From Table n we see 
that for electrons in the lens-shaped dot, the s-p energy 
level spacing Sgp — ep ~ Cs and p-d energy level spacing 
Spd = ~ are nearly equal, as assumed by the 2D 
harmonic model. The energy spacing Ssp and 5pd range 
from 50 - 80 meV (Fig. ^ , depending on the dot geome- 
tries. The electron energy spacings decrease with increas- 
ing QD base sizes. The electron energy spacing of alloy 
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TABLE I: Summary of the pseudopotential-calculated single-particle level spacing (in meV) of Ini_3;Gaa;As/GaAs quantum 
dots of different heights, base sizes, and Ga compositions, and hi are the i-th electron and hole single-particle energy levels. 
Ex is the lowest exciton energy. 
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Calculated energy-level diagram for 
lens-shaped InAs/GaAs quantum dot 
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FIG. 2: The schematic energy-level diagram (in meV) of 
a pure lens-shaped InAs/GaAs quantum dots, with height 
h=2.5 nm and base 6=20 nm. WLe and WLh denote the 
wetting layer energy levels for electrons and holes respectively. 
The CBM and VBM correspond to the conduction band min- 
ima and valence band maxima of (unstrained) bulk GaAs. Ex 
is the excitonic transition energy. 



dots are much smaller than those in pure InAs/GaAs 
QDs, because of reduced confinement. For Ga rich dots 
(x=0.3 - 0.5), the single-particle energy level spacings 
range from 30 - 45 meV. These values agree with the 
infrared absorption measurements^ of intra-band tran- 
sitions of alloy InGaAs QDs, which give Ssp ^41-45 
meV. When the Ga composition reaches x—0.5, the s- 
p energy level spacing 6sp becomes significantly different 
from the p-d energy level spacing Spd, thus deviating from 
harmonic potential approximation. 

(b) Shell definition: Figure |21 shows that the energy 
levels of electrons in a lens-shaped dot have well defined 
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(a) p-p splitting vs dot height 
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FIG. 3: Pseudopotential-calculated p level splitting of elec- 
trons and holes vs. (a) dot height and (b) dot base for 
InAs/GaAs dots. 



s, p d shell structure. However, while effective mass 
and k ■ p models predict degenerate p and d levels, for 
cylindrically-symmetric (e.g. lens-shaped) QDs, atom- 
istic calculations show that even in perfect lens-shaped 
dots, the p-p and d-d levels are split by 2 - 4 meV (Fig.|21) 
due to the actual C2v symmetry. We denote the two p 
levels as pi and p2, and similarly, the three d levels as c?i, 
^2 and c?3, in increasing order of energy. The results listed 



7 



and 5, 



dd 



^dn 



in Tabled show that 5pp = e^^ - 
are very sensitive to the aspect ratio of the dots while 
not being very sensitive to the alloy compositions. Fig- 
ure 121 depicts Spp — — vs. dot heights [Fig|3Ja)] 
and bases sizes [Fig|3Ib)]. In general, we see that 6pp in- 
creases with increasing dot height, and it decreases with 
increasing dot base size. 



2. Hole levels 

In contrast to electrons, hole single-particle levels (Ta- 
bic P) display a much more complicated behavior that is 
totally beyond the EMA description. 

(a) s-p spacing: As one can see from Tabled the hole 
s-p energy spacing ranges from 10 - 18 meV for dots of 
sizes we studied. These energy spacings are considerably 
smaller than those of electrons. The first confined hole 
state is found to be about 190 meV above the VBM of 
bulk GaAs, for the pure, h/b=2.5/20 nm dot (Fig. ^. 
Unlike the case for electrons, the energy spacing between 
hole s and p levels depends strongly on the height of the 
dots,^". while being relatively insensitive to the base size 
of the dots. 

(b) Shell definition: The well-defined s, p, d shell- 
structure for electrons does not exist for holes (Fig. 

as the pi-p2 and di-d2-d^ splitting are much larger than 
those for electrons The hole p-level splittings are also 
shown in Fig. O for different dot heights [Figl^Ja) ] and 
bases sizes [Figl^Jb)]. For the smallest dots, h/h — 2.5/25 
nm, the p splitting is about 11 meV (Table P), more than 
3 times the value for electrons. This splitting is about 
half of the hole s-p energy spacing. Note that the pseu- 
dopotential calculated p-p splitting is much lager than 1.3 
meV given by the k ■ p method (which includes piezoelec- 
tric effect)."'^ For taller dots, the p-p splittings are even 
larger. As a consequence, the p2 levels are energetically 
very close to the di levels, leading, as we will see below, 
to a nontrivial charging pattern that breaks Hund's rule 
and the Afubau principle.^^> 

(c) Wave function characters: An analysis of the wave 
function show that these levels have somewhat mixed 
S, P 01 D characters. For example, for a lens-shaped 
InAs/GaAs QD of height=3.5 nm and base=25 nm, the 
"s level" has 92% S character, and the two "p levels", 
Pi and p2, have 86% P character respectively (see Table 
II of Ref. l3lf) . Therefore, we can still label the single- 
particle levels as s, pi, p2, etc. These single-particle lev- 
els, s, pi and p2, do not have pure HH character either, 
being instead 91%, 86%, and 92% HH-like respectively. 
As the aspect ratio height/base increases, the mixture of 
angular momentum and HH and LH characters becomes 
stronger. For example, as shown in Ref. 43, for a In- 
GaAs/GaAs alloy dot, with 25.2 nm in base and 2 nm in 
height, s level has 90% S character, while both pi and 
P2 levels have 84% P character. When the height of the 
dot increases to 7.5 nm (with fixed base size), the leading 
angular momentum characters for these three levels are 



Coulomb energy of s orbitals 




1 2 3 4 5 6 7 
Dot height (nm) 

FIG. 4: Comparison of the electron-electron (ee) and hole- 
hole (hh) Coulomb energies of s orbitals Jas for lens-shaped 
dots with different heights. 



84% S, 78% P and 75 % P respectively. Similarly, for 
the flat dot (2 nm in height), the mixture of LH state 
character is about 4 - 9 %, but increases to 11 - 17 % for 
the tall dot (7.5 nm in height). 



B. Coulomb integrals: atomistic vs 2D-EMA 
description 

Another piece of information that decides the addi- 
tion energies [Eq. is the Coulomb integral between 
the particles [Eq. ifT^ ]. We list the pseudopotential- 
calculated Coulomb energies of s orbitals Jss for elec- 
trons in the first row of Table [H] and for holes in the 
first row of Table Iml For electrons, jiT^ ~ 22 - 25 
meV, and for holes jis '^ 20 - 27 meV for typical dots. 
These numbers can be directly compared with experi- 
mental value from electron/hole charging experiments, 
since Jss ~ A(lj2) [Eq. H12|) ]. The typical experimental 
values of Jss for electrons is about*' ^^'"''^ 19 - 27 meV and 
for holesiiiiS*^ 20 - 25 meV. We plot in FigH J^^ for elec- 
trons and holes vs. height for base 6 = 25 nm dots. For 
flat dots, the electron-electron Coulomb energy jis*^' is 
smaller than that of holes jis'*'. However, J^T^ is larger 
than jis'*' for taller dots. The crossover is at about 2.5 
nm for the 6 = 25 nm dots. Note that in our calculation, 
the two nearly degenerate electron p orbitals pi and p2, 
are spatially almost orthogonal to each other. However, 
in the simple 2D-EMA model,^^ the two degenerate p or- 
bitals p+ = Px + iPy and p_ — Px — iPy have same spatial 
function differing only by a phase factor. As a result, the 
exchange interaction between pi and p2 is much smaller 
than that of p+ and Furthermore, the simple 2D- 



EMA model predicts Jpj^p_ = Jp 
in the atomistic description, where Jj 



than Jp\p\ and Jp^p2 • 



PlP2 



which is not true 
is much smaller 



In a 2D-EMA model, there is only one free parameter 
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TABLE II: Electron- electron Coulomb energies Jap/ Jss and exchange energies Kap/Jss calculated numerically from atomistic 
pseudopotential theory for pure InAs/GaAs quantum dots compared with results obtained by the 2D-EMA model. The dots 
have base b= 25 nm with different heights. 
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TABLE III: Hole-hole Coulomb energies Ja^/Jss and exchange energies K^p/Jss calculated numerically from atomistic pseu- 
dopotential theory for pure InAs/GaAs quantum dots compared with results obtained by 2D-EMA model. The dots have base 
b— 25 nm with different heights. 
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in calculating the Coulomb integrals for each type of car- 
riers: the effective length for electrons and Ih for holes. 
This leads to simple geometric relationships as illustrated 
in Eq. (|17|l . These relations are listed in the second col- 
umn of Tables ITU and However, as noted in Sec. Ill Cl 
a detailed analysis of atomistic wave functiona^i*^ show 
that they are not of pure conduction band character for 
electrons or HH, LH characters for holes; nor do they 
have pure s, p angular momentum characters as predicted 
by 2D-EMA model. Since inter-valley, inter-band cou- 
pling and the underlying atomistic symmetries are also 
ignored in the 2D-EMA model, we might expect that 
the simple relations between the Coulomb integrals of 



Eq. I|17|) must be somehow different in the atomistic ap- 
proach relative to those predicted by the EMA, even for 
rather flat lens-shaped dots having parabolic-like energy 
level pattern, ep — ~ Cd — Cp, as chosen here. 

The ratios Jap/Jss and Kap/Jss derived from atom- 
istic calculations are compared with the 2D-EMA model 
relationships in Table|n]for electrons and in Table UTTl for 
holes for dots of different aspect ratios. To visualize the 
differences, we also plot the relative error for Coulomb 
energies ( Ja/3/ Jss)atm - {Jafs/ Jss)ema and for exchange 
energies (Kap/ Jss) c^tm - [Kais/ Jss)ema in Fig.|31(a) for 
electrons and in Fig. |S{b) for holes, where "atm" means 
atomistic results and "EMA" means the 2D-EMA model 
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FIG. 5: The relative error {Ja/s/ Jss)sutui — / Jss) jlma 

(closed circles) and {K^p / Jss)si\-m — (i^c/s/ Jss)ema (open cir- 
cles) for (a) electrons and (b) holes, vs. dot height for dots of 
base 6 = 25 nm. Here, "atm" means atomistic pseudopoten- 
tial calculation, and "EMA" means the 2D-EMA model. 



results. For flat lens-shaped dots which have a parabolic- 
like level spectrum, we see qualitative agreement between 
the atomistic calculations and the 2D-EMA model. For 
such flat dots, the errors are generally within 20% of Jss- 
For taller dots, the agreement is worse, as can be seen 
from Fig. [3 For even taller QDs (e.g., height/base > 
5/25 nm), the large biaxial strain will develop a "hole 
trap" at the InAs/GaAs interface, which lead to the hole 
localization on the interface of the dotsi^i*^ In these 
cases, where the hole (envelope) wave functions are to- 
tally different from those predicted by 2D-EMA model, 
the relationship between Coulomb and exchange energies 
are, accordingly, very different from those of 2D-EMA 
model. 

At first sight, it would seem that if the values of Jss 
in the 2D-EMA model are taken to be those calculated 
atomistically, the Coulomb integrals predicted by 2D- 
EMA model are only a few meV off from those calculated 
by atomistic theories. However, these small differences 
will change significantly the electronic phase-diagrams of 
carriers, as we will see below. 



C. Ground state configurations: atomistic vs. 
2D-EMA description 

1. Generic phase diagrams for 2D-EMA model 

To calculate the addition energies, we first need to 
know the spin and orbital configurations of the few- 
particle ground state, which minimizes the energy E{N) 
of the TV-particle system. For an illustration, we first use 
the "single configuration approximation" , i.e., retain only 
the Slater determinants that have the same orbital degree 
of freedom, and keep the exchange (spin-spin) coupling 
between these Slater determinants. The ground states 
are those configurations that have the lowest total en- 
ergy. For electrons, the calculated ground state for A^=l 
- 6 are shown in Fig. ^ which shows that the electron 
charging pattern follows Hund's rule. The ground state 
spin and orbital configurations for holes are listed in Ta- 
ble ^3 for N= 1-6 holes. In contrast to electrons, the 
hole charging patterns show complicated behaviors that 
defy the Hund's rule and the Aufbau's principle for most 
of the cases of N=5, 6 holes. 

To understand the differences between the ground- 
state configurations of electrons and holes and the driv- 
ing forces for these differences, we developed a gen- 
eral phase-diagram approachi^ that classifies the many- 
particle configurations for electrons and holes in quantum 
dots in terms of simple electronic and geometric parame- 
ters. To do so, we calculate for each particle number N, 
the configuration which minimizes the total-energy under 
the single-configuration approximation at different pi-p2 
splitting ((5pj_p2 = ~ ^pi) p2-di energy spacing 
(<5p2,di = Cdi — ep2) using Coulomb integrals J and K in 
units of Jss- This approach gives a phase diagram as a 
function of the parameters {N; Sp-^^p^, 
Jss, which yields for N=A two electronic phases 



P2i^P2,di} in units of 



'A = {sls^)iplpi) . 



(18) 



Here, we have adopted a spectroscopic notation for a sys- 
tem with cylindrical symmetry. For N — 5, there are 
three possible phases: 



= ish^){pl){pl){dl), 
= {sh-^)iplpi){pl) , 
= ish^)iplpi){dl) . 



(19) 



For N = 6, we find four phases, 

= is\s^)ipl)ipl)idl){dl 

3h = is\s^)iplp[)ipl)idl) 

= {sh^){plpi){plpi) , 

lE* = i.sh^){plpi){dldi) . 



(20) 



We first apply this approach to a 2D-EMA model. We 
relax the restrictions in 2D-EMA model of degenerate 
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FIG. 6: The most stable configurations in self-assembled InAs/GaAs dots, for N = 4, 5, 6 electrons/holes, as a function of the 
normalized pi-p2 splitting S{pi,p2) — ep^ — tp^ and p2-rfi splitting S{p2, di) ~ ^di — ep2- The notation of the configurations are 
given in Eqs. 1181 - 1201 . (a) Using the 2D-EMA model (the results apples to both electrons and holes) . (b) For electrons, using 
atomistic Coulomb and exchange integrals, (c) For holes, using atomistic Coulomb and exchange integrals. The 5{pi,p2) / Jss 
and S{p2,di)/Jss values of actual self-assembled dots are denoted by "e" and "h" in lenses for electrons and holes respectively. 
We also show the most stable electron configurations in the phase diagrams of a typical electrostatic dot represented by the 
circles for 2D-EMA model. 



shells (epi=ep2i ^rfi 



-Cp) and allow dp^^^p^ and Sp^ 



) and of equidistant shells 
to vary. The 



resulting phase diagrams are shown in Fig. IHIa) for iV= 
5 , 6. 



2. Ground states for specific dots in the 2D-EMA model 



To decide which of these phases is a ground state for 
a given dot, we need to know in Fig. EJa) the actual 
value of Spi^p2/Jss and Sp2^di/Jss for this realistic dot. 
For electrons in self-assembled dots, the single-particle 
energy spacing is usually more than twice the Coulomb 
energy,li& i.e., 5p2,di > '^Jss- For holes, (5p2,di = 
l.nJss was determined from recent experimentsiiii^ and 
5pi,p2 = is assumed in 2D-EMA model. This places 
in Fig. EJa), for both electrons and holes, phases '^S, 
^n, and as ground states for N—4, 5, 6, respectively, 
where the lens represents a lens-shaped self-assembled 
QD, and label "e", "h" inside the lens are symbols for 
electron and hole, respectively. As a comparison, we also 
show the electron ground-state phases for an electrostatic 
dot (~ 50 nm) represented by a circle. 



TABLE IV: Orbital occupations for holes in InAs/GaAs 
QDs of different height/base obtained from single-particle 
atomistic pseudopotential and many-particle CI calculations. 
The orbital occupations are given by the leading configura- 
tions in CI results. 
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3. Generic phase diagrams for atomistic approach 

We now use atomistic Coulomb and exchange energies 
J and K, listed in Table and IIIII to recalculate the 



phase diagram {N]5p^^pr^,5p^^di} for electrons and holes 
in Figs.inib) and Figs.EJc), respectively. We see that the 
phase boundaries change dramatically for both electrons 
and holes. For example, for 7V=6 electrons, phase ^S* 
disappears [FigEJb)], while phase disappear for N—h 
holes. 
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4- Ground states for specific dots in the atomistic approach 

We next use Sp-^^p^, Sp^.di and Jgs from Table HI and UTI 
to decide the ground states for electrons in self-assemble 
InAs/GaAs QDs. For example, for a InAs/GaAs dot 
with 2.5 nm in height and 20 nm in base, we have 
Jpj^pj/Jss ~ 0.1 and 6p2^di/Jss ~ 3. These parameters 
give the ground-state phase for 5 electrons, and the 
ground-state phase for 6 electrons. As we see, even 
though the phase boundaries calculated from atomistic 
theories are very different from those calculated from 2D- 
EMA model, the ground-state configurations are remain 
the same as those predicted by 2D-EMA model. This 
is because in the phase diagrams, the coordinates of the 
electrons in self-assembled dots are far away from other 
competing phases. For N=A electrons, the high spin state 
phase ■^E (Hund's rule) and low spin state are rela- 
tively close in the phase diagram. The intrinsic p-level 
splittings in the lens-shaped dots are about 1-4 meV 
(Table lU, which are smaller than the electron-electron 
exchange energies (~ 5 mV), and therefore the ground 
state is the high-spin phase ^E. However, if the dot 
shape is elongated, adding additional p-level splittings, 
the ground state could be low-spin phase ^A. Actually 
there are some experimental evidences showing that the 
ground state of 4 electrons could be a low spin stateAi^ 

In contrast to electrons, holes have Aeh < Jss, and 
large pi-p2 splitting, small p2-rfi splitting (Fig.O, which 
place the holes in a different region in the phase diagrams 
than electrons, where there are more competing phases. 
For example, for the dot with base b= 20 nm, and height 
h=2.5 nm, we have Sp-^^p^/Jss ~ 0.4 and Sp^^di/Jss ^ 
0.17. These relationships give, for this dot, the ground 
states '^E, ^A and ^E* for N=4, 5, and 6 holes, respec- 
tively, showing a nontrivial hole charging pattern that 
breaks the Aufbau principle.^® We also list the ground 
state configurations for dots of different geometries in 
Table II VI As we see, in most of the cases, the ground 
states are still ^E, ^A and ^E* as shown in Fig. El^c). 
However, there are some exceptions for very tall dots or 
very flat dots. For very flat dots (e.g., the h/b =2.5/27.5 
nm dot), the pi-p2 energy splitting (5.8 meV) is much 
smaller than the p2-di energy splitting (9.4 meV), These 
parameters place the ground state of N= 5 holes in phase 
^H and that of 6 holes in phase ^H of Fig.|S{c). For very 
tall dots (e.g., the h/b — 3.5/20 nm dot), the interfacial 
hole localization changes the hole-hole Coulomb integrals 
dramatically, and thus change the phase boundaries in 
[Fig. Etc)]. We find the ground states of the dot with 
h/b = 3.5/20 nm are ^H and ^E for 5 holes and 6 holes, 
respectively. 

5. Effects of configuration interaction 

If we use CI instead of single-configuration approxi- 
mation, the ground states are superpositions of different 
configurations, but the leading CI configurations have a 



significant weight, being 79%, 71% and 64% for 4, 5, 
and 6 holes, respectively, in the h/b = 2.5/20 nm dot. 
Since the weights of leading configurations are signifi- 
cantly larger than other configurations, we are justified 
in using leading configurations to represent graphically 
the ground states. It is worth noting that for very large 
electrostatically confined dots, where single-particle en- 
ergy spacing Se/ Jss 1 (down-left corner of the phase 
diagrams of Fig. ^ , the ground state mixes large num- 
ber of configurations, which have no significant leading 
configurations, and arc therefore in strongly correlated 
states*^ that are not discussed here. 



D. Calculated charging and addition energies 

Once we determined the ground state configurations, 
we can calculate the total energies using Eq. |0J. We 
calculate the ground state total energies for up to 6 elec- 
trons/holes for each dot. For electrons, in the CI ap- 
proach we used 6 single-particle electron levels (s, pi, p2, 
di, d2, dz) to construct all possible Slater determinants, 
while for holes, we used 8 single-particle hole levels. The 
total number of Slater determinants for 6 electrons is 924. 
For 6 holes, the total number of determinants is 8008. We 
plot in Fig. |7|the CI total energies for the ground state of 
6 carriers vs. number of single-particle states included in 
the CI expansions. The total energies converge to about 
1 meV if 6 single-particle states are used for electrons 
and 8 states are used for holes. 

The charging energies and addition energies are cal- 
culated using Eq. |^ and Eq. 10). The addition en- 
ergies calculated from CI calculations give results that 
are about 1 meV different from single-configuration re- 
sults for electrons, and 1-3 meV different for holes. 
The addition energies are summarized in Table Ex- 
perimentally, the charging energies of electronsii^ and 
holea^^i^ in self-assembled QDs are usually measured via 
capacitance spectroscopy, using gated structures.^* The 
(In,Ga)As/GaAs dots used in the electron charging ex- 
periments were roughly estimated to be 7 nm in height 
and 20 nm in base in Refs. m The measured A.(l, 2) ~ 
23 meV in Ref. ^ and - 21.5 meV in Ref. |1 which, 
as shown in Sec. ^ equals roughly JsT\ the electron- 
electron Coulomb interaction of s orbitals. The charging 
energies between s and o levels, Ae(2,3) was estimated 
to be 57.0 meV in Ref. |^ The average addition energies 
between p-states are estimated to be 18 meV in Ref. |^ 
and 14 meV in Ref. Furthermore, in Ref. 0, Ae(4, 5) 
is almost twice as large as Ae(3,4) and Ae(5,6), which 
might be related to the breakdown of Hund's rule as a 
consequence of irregular shape of the dots. These results 
are listed in Table and compared with our theoreti- 
cal results. We see that the electron addition energies 
of an (In,Ga)As/GaAs dot, with b/h = 3.5/25 nm and 
Ga composition x— 0.15, agree very well with the above 
experimental results, which show Ae(l,2)= 20.5 meV, 
Ae(2,3)= 61.2 meV, and the average addition energies 
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TABLE V: Summary of atomistically calculated electron addition energies Ae and hole addition energies (in meV) of 
Ini-iiGa^ As/GaAs quantum dots for various dot heights, base sizes and Ga compositions. The experimental addition energies 
for electrons are extracted from Ref. iQ, and the addition energies for holes are taken from Ref. ITU 
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FIG. 7: CI total-energy convergency test for (a) six-electron 
and (b) six-hole states with increasing number of single- 
particle orbitals in the bases. The test is done for a pure 
InAs/GaAs dot, with base 6=25 nm, and height h—3.5 nm. 



between p-states of about 16 meV. 

The experimental hole addition energies are taken from 
Ref. 0, which gives Ah(l,2)=23.9 meV, comparable 
to that of Ae(l,2). However, the addition energy be- 
tween s and p orbitals, A/i(2, 3)=34.2 meV, is signifi- 
cantly smaller than Ae(2,3) ~ 57 meV. This result re- 
flects that the s-p energy spacing of holes is much smaller 
than that of electrons. As seen from Table IVl our calcu- 
lated addition energies of pure and flat (height=2.5 nm) 



InAs/GaAs dots agree very well with this experiment. 

To study trends of addition energies for the electrons 
and the holes, we depict the electron addition energies 
for different dot heights [Fig.|Sfa)], bases [Fig.|SIb)] and 
alloy compositions [Fig. ISfc)]. Similarly, the hole addi- 
tion energies are plotted in Figs.|Hl(a),|H|[b), andl^c) for 
dot heights, bases and alloy compositions. We see the 
following: 

(i) Up to six carriers, the electron addition energies 
A(A^-1, TV) have a single peak at (A^-l)=2. The peak is 
due to the single-particle energy gap between the s-shell 
and p-shell. On the other hand, all hole addition energies 
have two peaks at {N — 1)=2 and {N — 1)=4 respectively, 
where the first peak come from the single-particle energy 
gap between the s orbital and pi orbital and the second 
peak is associated with the energy difference between the 
Pi and higher energy orbitals. 

(ii) Electron addition energies decrease with increasing 
height and base of the dots, and the hole addition ener- 
gies share the same trend. However, the electron addition 
energies are more sensitive to the base of the dots and 
relatively insensitive to the height. In contrast, the hole 
addition energies are very sensitive to the heights of the 
dot, and relatively insensitive to the base. This depen- 
dence suggests that the electron wave functions are more 
2-dimensional like in the QDs, while hole wave functions 
are more extended in [001] direction. 

(iii) Electron addition energies show a simple trend for 
alloyed Ini„2,Gaa;As/GaAs QDs, which are consistently 
smaller for Ga rich dots. However, the hole addition ener- 
gies of alloy dots show more complicated behaviors. The 
reasons of this complication are the following: (1) alloy 
dots have different trends in hole single-particle energy 
level spacings (see Tabled; (2) the ground state configu- 
rations might be different for different alloy compositions. 
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FIG. 8: Electron addition energies of (a) pure InAs/GaAs 
dots with base 6=25 nni and height h=2.5 nni and 3.5 nrn; 
(b) pure InAs/GaAs dots with height h=3.5 nni and base 
6=20 nm, 25 nm and 27.5 nm; and (c) Ini-^Ga^As/GaAs 
dots with height h=3.5 nm and base 6=25 nm at x=0, 0.15, 
0.3 and 0.5. 



FIG. 9: Hole addition energies of (a) pure InAs/GaAs dots 
with base 6=25 nm and height h—2.5 nrn and 3.5 nrn; (b) pure 
InAs/GaAs dots with height h=3.5 nm and base 6=20 nm, 25 
nm and 27.5 nm; and (c) Ini-xGa^As/GaAs dots with height 
h=3.5 nm and base 6=25 nm at x=0, 0.15, 0.3 and 0.5. 



IV. SUMMARY 

We systematically studied the electron/hole addition 
energy spectra using single-particle pseudopotcntial plus 
many-particle CI methods. Considering the single- 
particle step, we find for electrons that there is a shell 
structure and that the p-p and d-d splittings are about 
1-4 mcV depending on the dot geometry. For holes, 
the single-particle step reveals large (5 - 11 meV) pi- 
P2 splitting and absence of a well defined shell struc- 
ture. Considering the e-e and h-h Coulomb integrals, we 
find that atomistically calculated ratios between various 
Coulomb integrals Jij differ by about 20% from those in 
the 2D-EMA model calculation. These differences lead 



to many-particle phase diagrams that differ significantly 

from those predicted by the 2D-EMA model. In par- 
ticular, the "unusual" hole single-particle spectrum and 
Coulomb integrals lead to many-particle ground states 
that defy the Hund's rule and the Aufubau principle for 
holes. The predicted ground-state configurations and the 
addition energies calculated in this pseudopotential plus 
CI scheme compare well with experiments. 
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